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Abstract 



We perform a detailed study of the decays of the heavier r slepton (f2) and r- 
sneutrino (t'r) in the Minimal Supersymmetric Standard Model (MSSM). We show 
that the decays into Higgs or gauge bosons, i.e. ff + or 

— > z^T- + (H~ or W~), and t'r + {H~^ or W^"*"), can be very important due 

to the sizable r Yukawa coupling and large mixing parameters of f . Compared to 
the decays into fermions, such as f2^T~ + 'j^ and f2 i^t + Xj ^ these bosonic 
decay modes can have significantly different decay distributions. This could have 
an important influence on the search for f2 and Ur and the determination of the 
MSSM parameters at future colliders. 



In the Minimal Supersymmetric Standard Model (MSSM) |^ supersymmetric (SUSY) 
partners of all Standard Model (SM) particles with masses less than 0(1 TeV) are in- 
troduced. This solves the problems of hierarchy, fine-tuning and naturalness of the SM. 
Hence discovery of all SUSY partners and study of their properties are essential for test- 
ing the MSSM. Future colliders, such as the Large Hadron Collider (LHC), the upgraded 
Tevatron, e~^e~ linear colliders, and colliders will extend the discovery potential for 

SUSY particles to the TeV mass range and allow for a precise determination of the SUSY 
parameters. 



In this article we focus on the sleptons of the third generation, i. e. staus (fi 2; < rrif 
and tau-sneutrino (z/t-). These particles may have properties different from the sleptons of 
the other two generations due to the sizable r Yukawa coupling. Production and decays 
of fi and z/t- were studied in [|, ^, ^. Like other sleptons, they can decay into fermions, 
i. e. a lepton plus a neutrahno (x^) or chargino (xf)'- 

Xj , t^T-^r-xt , 

with i,j = 1,2 and k = 1,...,4. In addition, the heavier stau f2 and the tau-sneutrino 
can also decay into bosons ||, ^, i. e. a lighter slepton plus a gauge boson 



(2) 



or a Higgs boson 



- ff (/^o, ffo. A') , 

T2 ^ Ur H , i/^ Ti . 

The decays in Eqs. (|^) and (0) are possible in case the mass splitting between the sleptons 
is sufficiently large. 

In the present article we perform a more general analysis than [0, We point out that 
the f2~ and Ur decays into gauge or Higgs bosons of Eqs. (H) and (H) can be very important 
in a large region of the MSSM parameter space due to the sizable r Yukawa coupling and 
large f mixing parameters. This importance of the bosonic modes relative to the fermionic 
modes of Eq. (|^) could have a significant infiuence on searches for fg" and 9^ at future 
colliders. An analogous study for the heavier stop (^2) and sbottom (62) was performed 
in §,§. 

First we summarize the MSSM parameters in our analysis. In the MSSM the stau 
sector is specified by the mass matrix in the basis (f^, f^) |0, ||] 

= ( < ^^T^ ] (4) 

with 



m?^ = M| + m| cos 2/3 (sin^ Ow - h + ml, (5) 
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m^^ = Ml; -ml cos 2/5 sin^ 9w + m^, (6) 
arm^r = (v4t- — /itan/5) mr . (7) 

Ml ^ and A^- are soft SUSY-breaking parameters, fi is the higgsino mass parameter, and 
tan/3 = V2/V1 with t>i (t>2) being the vacuum expectation value of the Higgs field {H2). 
Diagonalizing the matrix (^) one gets the mass eigenstates fi = f2,cos6'f + fRsin6'f, and 
f2 = — fi:,sin6'f + TRCOsOf with the masses mf^, ("^n < J^fa) ^^^1 the mixing angle 
Of- The stau mixing is large if |m|^ — m?^| < |aT-mT-|, which may be the case for large 
tan /3 and /i. The mass of Vr is given by 

m%=Ml + Knlcos2[3. (8) 

The properties of the charginos xf {i = 1,2; m-± < ""^^i) and neutralinos Xfc = 
1,...,4; m^o < ... < TM^o) are determined by the parameters M, M', /i and tan/3, 
where M and M' are the SU(2) and U(l) gaugino masses, respectively. Assuming gaug- 
ino mass unification we take M' = (5/3) tan^ The masses and couplings of the 

Higgs bosons H^, A^, and if^, including leading radiative corrections, are fixed by 
rriA, tan/3, fi, mt, mb, Mq, Mjy, Mfj, At, and Af,, where Mq jy and At^f, are soft 
SUSY-breaking parameters in the (t, b) sector. (h^) and A^ are the heavier (lighter) 
CP-even and CP-odd neutral Higgs bosons, respectively. For the radiative corrections to 
the and masses and their mixing angle a we use the formulae of Ref . [^] ; for those 
to mH+ we follow Ref. Il0|.[| We treat M^ij^^ gjjjjy and ^{r,t,f)} as free parameters. 



The widths of the ii = T2 or Ur decays into Higgs and gauge bosons are given by {k 
l,-,4) i: 



m - H,) = {G,,u)\ r(£, ty V) = 3 (Q,y)'. (9) 



Here = or Vr (with the indices i and j to be omitted for Vr), = {h^, H^, A^, H^} 
and V = {Z^,]^^}. k^^x = '^("^f , ^x) is the usual kinematic factor, K{x,y,z) = 

J 

(x^ + y"^ + z'^ — 2xy — 2xz — 2yzY^'^. Notice an extra factor k? /m\^ for the gauge boson 
modes. Gijk denote the slepton couplings to Higgs bosons and Cijv those to gauge bosons. 
Their complete expressions, as well as the widths of the fermionic decays, are given in |Q . 

Since is rather small, we need a large difference between m^^ ~ m^^ and m^^ in order 
to realize a mass splitting between (fi, f2, Vr) large enough to allow the bosonic decays 
in Eqs. (2, 3). In this case the fi — mixing is rather small. Thus, in this article 
we consider two patterns of the mass spectrum of the sleptons: m^^ < mf^ ~ m^^ with 
(ri,f2) {fR,fL) for > mrj,, and m^^ > m^^ ~ m^^ with (fi,f2) ~ {tl^tr) for 

^ Notice that p, nO| have a sign convention for the parameter /i opposite to the one used here. 
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rufj^ < ruf^. The bosonic decays are therefore basically the decays of {jl^Vt) into tr or 
vice versa. This is in strong contrast to the case of the t2 and 62 decays 0. In the latter 
case, large mass splittings can be obtained even for Mq ~ ~ due to the large 
top and/or bottom Yukawa couplings, which can also cause large left-right mixings and 
a complex decay spectrum. 

The leading terms of Gijk and Cijv which are relevant for the bosonic decays are given in 
Table 1. Here the Yukawa coupling hj. is given by 

hr = g rrir / {V2 mw cos (3) . (10) 

Higgs bosons couple mainly to {jl,Ut) — fji combinations. These couplings are pro- 





~ /It- (/icosa + sina) cos26'f 




~ h^- {fi sin a — At- cos a) cos 26 f 




~ hr {fi cos P + At- sin P) 




~ g sin26'f 




~ hr (/i cos P + Ar sin P) sin Of 




~ hj. (/i cos P + At- sin P) cos Of 




~ g cos^f 




~ —g sin6'f 



Table 1: Leading terms of the slepton couplings to Higgs and gauge bosons. 

portional to the Yukawa coupling /i^ and the parameters Aj. and yU, as can be seen in 
Table 1. Notice the factor cos26'f in the fj^fg/i" and f^f2H^ couplings, but not in the 
fj^fgA'' coupling. In our case | cos2^^f | ~ 1 due to the small — tr mixing, unlike the case 
of ^2 and 62 decays 0. Similarly, the VrfiH^ couphng is not suppressed by f mixing for 
~ tr. Hence the widths of the decays into Higgs bosons can be large for large tan/?, 
Ar, and /x. Notice that the — mixing enhances the mass splitting between fi and f2, 
which results in a larger phase space for the bosonic decays of f2. In contrast, the gauge 
couplings which are relevant for the bosonic decays are suppressed by the small mixing 
of sleptons, since the gauge interactions preserve the chirality of sleptons. However, this 
suppression is largely compensated by the extra factor /nty in Eq. (^). In fact, since 
the gauge bosons in the decays (|[) are longitudinally polarized, the widths of the decays 
into gauge bosons for m^. — m^(i) ^ my are approximated by those into the corresponding 

Nambu-Goldstone (NG) bosons. As a result, the decays into gauge bosons are enhanced 
when the couplings of the sleptons to the NG bosons (oc arrrir/mv) are large. On the 
other hand, the fermionic decays are not enhanced for large [A^, /i, tan/5). Although the 
slepton-Higgsino couplings increase with tan/?, their contribution to the kinematically 
accessible decays (|1|) is suppressed for > rrif^. Therefore the branching ratios of the 
bosonic decays (|^, |) are expected to be large for large A^, /x, and tan/3. 
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We now turn to the numerical analysis of the f2 and Ur decay branching ratios. We 
calculate the widths of all possibly important two-body decay modes of Eqs. (|I|, ^ 
Three-body decays are negligible in this study. We take m,- = 1.78 GeV, rrit = 175 GeV, 
rrib = 5 GeV, = 91.2 GeV, sin^ 6^^^ = 0.23, m^r = rrizCosOw and a;(m^) = 1/129. 
In order not to vary too many parameters we fix M = 300 GeV, uia = 150 GeV, and 
Mq = Mfj = Mjy = At = Af) = 500 GeV for simplicity. In our numerical study we 
take m^j, rrif^, A^-, /i, and tan/? as input parameters. Note that for a given set of the 
input parameters we have two solutions for (M^, M^) corresponding to the two cases 
rufj^ > rrif^ and m^;^ < rrifj^. In the plots we impose the following conditions in order to 
respect experimental and theoretical constraints: 



(i) m-± > 100 GeV, m,>^ > 45 GeV, rriho > 90 GeV, m-^ i^ i^ > m^o, m^o > 80 GeV, 



ii) ^2 <3(M| + M| + m^J, < 3(M?+M?+m2,J, and < 3 (M? +M|, + m2,J, 
(m 



where m|^^ 



\ + m?z) sin^ P — \ ^| and m^^ 



-|- m|) cos^ /? — I m 



(iii) Ap{t-b) < 0.0012 Q using the formula of p 



Condition (i) is imposed to satisfy the experimental bounds on xt, Xi; ^i b, g, and 
from LEP and Tevatron \T^\. Note that m^o > 80 GeV is imposed in order to evade 



the experimental bounds on m^^ and . Condition (ii) is approximately necessary to 
avoid color and charge breaking global minima ||15| and to exclude unrealistically large At- 
Condition (iii) constrains fi and tan (3 in the squark sector. We note that the experimental 
data for the b s'j decay give rather strong constraints fl^ on the SUSY and Higgs 
parameters within the minimal supergravity model, especially for large tan (3. However, 
we do not impose this constraint since it strongly depends on the detailed properties of 
the squarks, including the generation-mixing. 

In Fig. 1 we plot in the A^-fi plane the contours of the branching ratios of the Higgs 
boson modes BR(^ ^ i'H) = ^ BR[ £ ^ f + {h°, ^4°, H^)], the gauge boson modes 
BR{i i'V) = j:BR[£^i' + {Z^, W^) ], and the total bosonic modes BR{£ £'B) = 
BR(£ ^ i'H) + BR{i i'V), with i = {f^, z>,) and i' = (ff , z>,). We take m^, = 250 
GeV, rrif^ = 500 GeV, tan/3 = 30 and show two cases m^^ < rrif^ and m^^ > rrif^. Note 
that i>r decays into bosons and decays into Ur are kinematically forbidden if rrif^ < mf^ 
and > m^^, respectively. 

We observe that BR{i i'H) increases with l^,-! while BR{i i'V) increases with 
\ix\. This dependence on (A^, /x) is explained as follows. Note first that for large tan/? 
the mixing between Hi and H2 is rather small. Hence, for ttla > friz, {H^, A^, H^) are 
mainly Hi while h'^ and the NG bosons are mainly i^2- The couplings of (f, i/^) to Hi 
and H2 are ~ hrAr and hrfi, respectively. Therefore the decays to {H^, A^, H"^) are 
enhanced for large l^^l, whereas those to and to gauge bosons are enhanced for large 
1/^1 . This property can also be derived directly from Table 1 by noting that | sina| <C 1 
and sin/3 ~ 1. As a result, the total bosonic branching ratio BR{i — > i'B) becomes large. 
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and even dominant for mf^ < mr^, in a wide region of the Ar^n plane, especially for large 
\At-\ and/or as seen in Fig. 1. We also see that the branching ratios of the bosonic 
decays are almost unchanged by A^. —A^ and/or /i — /x. 

In Fig. 2 we show the individual branching ratios of the f2 and Vr decays as a function of 
tan 13 for m^, = 250 GeV, nif^ = 500 GeV, A^ = 800 GeV and /i = 1000 GeV, for the two 
cases > and m^^ < nir^. We see that the branching ratios of the boson modes 
increase with tan/3 and become important, and even dominant if rrif^ < m^^, for large 
tan/3(>15). As already explained, this comes from the increase of h^- and 0^777.7- with 
tan/5. We have checked that r(f2 tiH^) itself increases with tan/3. 

In Fig. 3 we show the dependence of the 7^2 and 9^ decay branching ratios for rrif^ = 250 
GeV, Ar = 800 GeV, = 1000 GeV and tan /3 = 30. We see that the branching ratios 
of the bosonic decays decrease with increasing rrir^. This behavior comes from the fact 
that in the large nif^ limit the decay widths of the bosonic and fermionic modes are 
proportional to rrif^ and nif^, respectively. 

In all figures one can see that the total branching ratio of the bosonic decays of f2 is 
substantially larger for ruf^ < nif^ than for nif^ > rrif^. The reason is the following: 
First, as already explained, the decays of f2~ into Ur are kinematically allowed only for 
?7if^ < rufj^, as seen in Figs. 2 and 3. Second, for the parameter set used in Figs. 2 and 
3, Xi ~ W~ . Hence the decay f2 — > J^tXi is strongly suppressed for nif^ < rrif^ (i.e. 
T2 ~ tr), but not for m^^ > m^^, as seen in Figs. 2 and 3; the f2 decay into x-i ~ H~ 
is kinematically forbidden. This results in a rather strong enhancement of the bosonic 
modes for the case m^^ < m^^. 

We find that the importance of the bosonic modes is fairly insensitive to the choice of 
the values of rriA, M, Mqjj f) and At^t- The decays to H^, A'^ and are kinematically 
suppressed for large ttia- However, the remaining and gauge boson modes can still 
be important and even dominant. For example, for m^^ = 250 GeV, rrif^ = 500 GeV, 
A^ = 800 GeV, /i = 1000 GeV, tan/3 = 30 and ruf^ < rrif^, the fa decays into A°, 
H^) are forbidden for ttia > 260 GeV. Nevertheless, in this case we have BR(f2 fih^) 
= 13% and BR(f2 — > iV) = 52%. We have also checked that our results do not change 
significantly for smaller values of M in the range where the condition (i) is satisfied. 

Now we discuss the signatures of the fa and Ur decays. We compare the signals of the 
decays into bosons (Eqs. |^)) with those of the decays into fermions (Eq. (|l|)). The 
bosonic decays always produce cascade decays. In principle, the final states of the bosonic 
decays can also be generated from fermionic decays. For example, the final particles of 
the decay chain 

^ fr + ih', H', A' or Z°) ^ (r-x?) + (bb) (11) 

are the same as those of 

^2 ^r-+ xlsA ^r- + ((/7°, H', A' or + x?) - + (bbxi)- (12) 
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Likewise, 

Vr + {H^ 01 W+) (r-Xi) + iqq') (13) 

has the same final particles as 

i>r^r- + xt2 ^r- + or W+) + x?) ^ + {qq'xD- (14) 

Nevertheless, the decay distributions of the two processes (|1^) and (p!^ ) (([I3|) and (p!4D) are 
in general different from each other due to the different intermediate states. For example, 
the T~ in the chains ([TT]) and ( [T^ ) tends to be softer than the r~ in (|T2|) and (|n|), respec- 
tively. A similar argument holds for the quark pairs in the decay chains. Moreover, the 
distribution of the missing energy-momentum carried by x? could be significantly different 
in ( pTD and (|12[) ((p!3D and ([l^)) since it is emitted from a different sparticle. Detailed 
Monte Carlo simulations are necessary to investigate the experimental consequences of 
the bosonic decays. 

In conclusion, we have shown that the decays of f2 and v-r into Higgs or gauge bosons, 
such as f2 —>■ fi + (/i°, or Z^) and fg' ^ i/^ + {H^ or W~)^ can be very important 

in a fairly wide MSSM parameter region with large mass splitting between fi and f2, 
large tan/3, and large 1^4^-1 and/or Compared to the fermionic decay modes, these 
bosonic decay modes could have significantly different decay distributions. This could 
have important implications for the searches of f2 and and the determination of the 
MSSM parameters at future colliders. 
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Figure Captions 



Figure 1: Branching ratios of f2 and Ur decays in the A^-n plane for = 250 GeV, 
= 500 GeV and tan /3 = 30 in the cases of rrif^ < mf^ (a - c) and mf^ > mf^ (d - i) 
; (a, d) EBR[f2- ^fr + (/?o,i/o,AO), Ur + H-], (b, e) EBR[f^ ^ fr + Z*", u^ + W-], 
(c, f) EBR[f2- ^rr + (/iO,ifO,AO,Z"), i), + iH-,W-)], (g) BR[ ^ ff + /J+ ], (h) 
BR[ rf + W^], and (i) E BR[ t/^ ^ rf + 1^+) ]. The gray areas are excluded 
by the conditions (i) to (iii) given in the text. 

Figure 2: tan (3 dependence of f2 (a, b) and Ur (c) decay branching ratios for m^^ = 250 
GeV, = 500 GeV, Ar = 800 GeV and = 1000 GeV in the cases of rrifj^ < nif^ (a) 
and > m^^ (b, c). "Gauge/Higgs + X" refers to the sum of the gauge and Higgs 
boson modes. The gray areas are excluded by the conditions (i) to (iii) given in the text. 

Figure 3: nif^ dependence of f2 (a, b) and Dj. (c) decay branching ratios for ruf^ — 250 
GeV, At — 800 GeV, ji — 1000 GeV, and tan /3 = 30 in the cases of vrif^ < m^^ (a) and 
ruf^ > ruf^ (b, c). "Gauge/Higgs + X" refers to the sum of the gauge and Higgs boson 
modes. The gray areas are excluded by the conditions (i) to (iii) given in the text. Note 
that mo^ ~ in (c). 
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Fig.l 
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Fig.2 
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